Cold atoms in non-Abelian gauge potentials: 
From the Hofstadter "moth" to lattice gauge theory 
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We demonstrate how to create artificial external non-Abelian gauge potentials acting on cold atoms in optical 
lattices. The method employs n internal states of atoms and laser assisted state sensitive tunneling. Thus, 
dynamics are communicated by unitary n x n-matrices. By experimental control of the tunneling parameters, the 
system can be made truly non-Abelian. We show that single particle dynamics in the case of intense U (2) vector 
potentials lead to a generalized Hofstadter butterfly spectrum which shows a complex "mofh"-like structure. We 
discuss the possibility to employ non-Abelian interferometry (Aharonov-Bohm effect) and address methods to 
realize matter dynamics in specific classes of lattice gauge fields. 



One of the most significant trends in the physics of 
ultra-cold gases nowadays concerns, without any doubts, 
strongly correlated systems. Apart from systems at BCS- 
BEC crossover, low dimensional systems, and atomic lattice 
gases, perhaps the most fascinating possibilities are offered 
by atomic gases subject to the effects of artificial magnetic 
fields. A feasible realization, albeit experimentally challeng- 
ing, is an atomic system in a rotating trap. Trap rotation mim- 
ics the effects of external homogeneous magnetic fields. In 
the situation, when the rotation frequency approaches the trap 
frequency, the systems may allow to study the fractional quan- 
tum Hall effect (FQHE) and Laughlin liquids 1 1]. Several ex- 
periments with gases in rotating traps are currently underway 
1211 . In this context, rotating dipolar lattice gases (for a review, 
see 1 3]) offer new possibilities and advantages compared to 
gases with short ranged interactions, as the energy gap be- 
tween the Laughlin state and quasi-hole excited states survives 
the large TV limit. This allows for the realization of Laughlin 
liquids with a few hundred atoms at realistic temperatures |4]. 
Very recently, the creation of artificial magnetic fields in terms 
of electro-magnetically induced transparency (EIT) has been 
proposed in Ref. |8]. 

Alternatively, magnetic field effects can be realized in a 
lattice by introducing appropriate phase factors for tunnel- 
ing amplitudes. If multiplied around an elementary plaquette, 
these lead to phases proportional to the magnetic flux pen- 
etrating the plaquette. Jaksch and Zoller, and subsequently 
other groups, have recently proposed methods to realize such 
"artificial" magnetic field effects in lattice gases. These em- 
ploy atoms with multiple internal states, laser assisted tunnel- 
ing, lattice tilting (acceleration), and other experimentally ac- 
cessible techniques 0, S 01- The physical features in such 
"artificial" magnetic fields are extremely rich. For single 
atoms, the spectrum exhibits the fractal Hofstadter "butter- 
fly" structured] illustrated in Fig. In weakly interacting 
or weakly disordered systems, the modifications of the but- 
terfly due to interactions and disorder, respectively, may be 
studied |5]. Finally, considering strong magnetic fields in the 



FIG. 1 : The Hofstadter butterfly. Eigenenergies in natural units plot- 
ted versus magnetic flux per plaquette a — p/q,£ [0,1], with 
q < 100. The lowest Bloch band splits into subbands and forms 
a fractal for the set of all q's. 
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limit of strong interactions, Laughlin like states are expected 
0|. In this Letter, we show that the use of atoms with more in- 
ternal states, the application of state dependent, laser assisted 
tunelling (cf. IIOII ). and coherent transfer between internal 
states, allows for a generalization of the above methods and 
creation of "artificial external magnetic fields" corresponding 
to non-Abelian U (n), or SU(n) gauge fields 1 1 1 ]. In this case, 
tunneling amplitudes are replaced by unitary matrices whose 
product around a plaquette is non-trivial and its mean trace 
(Wilson loop) is not equal to n To our knowledge, the 
physics of matter in designed non-Abelian gauge fields has not 
been studied at all in the AMO physics context. Although sin- 
gle particles in SU (n) monopole fields have been intensively 
studied in high energy physics fTjll . other field configurations 
have not attracted such interest. Specific examples of external 
gauge fields have also been discussed with respect to NMR, 
nuclear and molecular physics fl4ll . Recently, there has been 
a proposal to create non-Abelian gauge fields using EIT fl5ll . 
However, all of these proposals either do not deal with inter- 
acting many particle systems in such gauge fields, or they do 
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not consider lattice gauge fields of ultra-high strength. 

Having introduced the designed external non-Abelian 
gauge fields in optical lattices, we study their influence on sin- 
gle particle dynamics. We generalize the Hofstadter butterfly 
to the non-Abelian case and show that in case of intense U (2) 
gauge fields, the spectrum (as a function of two "magnetic 
fluxes") exhibits a complex "moth" structure of holes. We 
discuss the possibility to observe the U (2) Aharonov-Bohm 
effect, to realize non-Abelian interferometry, and to simulate 
lattice gauge dynamics using random field configurations. 

We start from the physics of a single particle in a two- 
dimensional square lattice of spacing a in the presence of an 
Abelian magnetic field B. When the lattice potential is suffi- 
ciently strong, the single particle non-relativistic Hamiltonian 
in the tight binding approximation is given by the so-called 
Peierl's substitution [9], and reads 
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where Vq is the optical potential strength, p is the momen- 
tum operator, and A is the magnetic vector potential. Given 
[pj ,Aj] — 0, the single-particle wave equation reads 
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+e A »ip(x,y+a) + e 1 ?^ y ip(x,y—a) = — — <ip(x, y) 

Vq 

The choice of A determines the magnetic field B, which in 
turn determines the behavior of the system. For B — Be z , one 
may choose the potential A = (0, Bx, 0), thus, solely tunnel- 
ings in the y-direction acquire phases. This makes the prob- 
lem effectively one-dimensional and (|2) can be transformed 
into Harper's equation I17II : 

g(m+l) — g(m—l) + 2cos(27rma — v)g(m) = eg(m) (3) 

by using the ansatz ip(ma, no) — e lvn g{m), where x = ma, 
y = na, and e = —E/Vq. Solving the eigenvalue problem for 

2 p 

specific, i.e., rational values of the magnetic flux a = °° c 
per elementary plaquette makes the problem periodic. This 
results in a band spectrum whose gaps form the famous Hofs- 
tadter butterfly (Fig.Q. Note that the regime of this spectrum 
requires finite values of a, i.e., magnetic fields B ~ 1/a 2 , 
which in the continuum limit a — > become ultra-intense. A 
physical system satisfying Q can be realized in a simpler way 
in a ID optical lattice setup with standard tunneling and an 
overlaid super-lattice potential with spatial period 1/a. This 
statement also holds for the non-Abelian generalizations in 
the scope of this manuscript. Anyhow, the full 2D realization 
proposed in [5] allows for a more powerful control opening 
up completely novel features in such a system. 

We consider an atomic gas in a 3D optical lattice and 
assume that tunneling is completely suppressed in the in- 
direction, so that, effectively, we deal with an array of 2D 
lattice gases and are able to restrict ourselves to one copy. 
The atoms occupy two Internal hyperfine states \g), |e), and 



the optical potential traps them in the state \g) and |e), re- 
spectively, in every second column, i.e., for the y coordinate 
equal to . . . , n — 1, n + 1, . . . (. . . , n, n + 2, The result- 
ing 2D-lattice has thus the spacing A/2 (A/4) in the x- (y-) 
direction. The tunnel rates in the x direction are due to ki- 
netic energy; they are spatially homogeneous and assumed 
to be equal for both hyperfine states. The lattice is tilted in 
the y-direction, which introduces an energy shift A between 
neighboring columns. Tilting can be achieved by accelerat- 
ing the lattice, or by placing it in a static electric field linearly 
dependent on y. 

By this, standard tunneling rates due to kinetic energy are 
suppressed in the y direction. Instead, tunneling is laser as- 
sisted, and driven by two pairs of lasers resonant for Raman 
transitions between \g) and |e), i.e., n <-> n ± 1. This can 
be achieved because the offset energy for both transitions is 
different and equals ±A. Detunings of the lasers are chosen 
in such a way that the effect of tilting is cancelled in the ro- 
tating frame of reference. The lasers consist of running waves 
in the ±x-direction, so that the corresponding tunneling rates 
acquire local phases cx-p(±iqx). 

In order to realize "artificial" non-Abelian fields in a similar 
scheme, one has to use atoms with internally degenerate Zee- 
man sublevels in the hyperfine ground state manifolds, \gi), 
and |e,) with i = 1, . . . , n, whose degeneracy is lifted in ex- 
ternal magnetic fields. These internal states may be thought 
of as "colors" of the gauge fields. Promising fermionic can- 
didates with these properties are heavy Alkali atoms, for in- 
stance, 40 K atoms in states F = 9/2, M F = 9/2, 7/2, . . ., 
and F = 7/2, M F = -7/2, -5/2, . . .; in particular, they al- 
low for realizing "s pin" dependent lattice potentials and spin 
dependent hopping 1 1^1 • 

Having identified the "colors", one should essentially re- 
peat the scheme of Ref. @] with two additional modifica- 
tions: laser assisted tunneling rates along the y-axis should 
depend on the internal state, although not necessarily in the 
sense of Ref. flflll . They rather should be performed in such a 
way that for a given link \gi) to |e»), tunneling should be de- 
scribed by a non-trivial unitary tunneling matrix U y (x) being 
a member of the "color" group (Z7(n), SU(n), GL(n) etc.). 
For unitary groups, the tunneling matrix U y (x) can be rep- 
resented as cxp(iaA y (x)) . Here, a is real, and A y (x) is a 
Hermitian matrix from the corresponding gauge algebra, e.g., 
u(n) or su(n). Since transitions from \g.i) to |ej) correspond 
to different frequencies for each i, they are driven by differ- 
ent running wave lasers. Thus, they may attain different phase 
factors exp(±iqix). 

In order to create gauge potentials that cannot simply be re- 
duced to two independent Abelian components, tunneling in 
the x-direction should also be Raman laser assisted and should 
allow for coherent transfer between internal Zeeman states. 
We will assume for now that this is proceeded by the same uni- 
tary tunneling matrix U x for both hyperfine state manifolds, 
although more general situations are feasible and basically in- 
teresting. To assure a genuine non-Abelian character of the 
fields, it is necessary that [U x , U y (x)] ^ 0. We stress that all 
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FIG. 2: Optical lattice setup for U(2) gauge fields: Red and blue 
open semi-circles (closed semi-circles) denote atoms in states \gi) 
and 1 32), respectively (|ei) and 62)). Top) Hopping in the re- 
direction is laser assisted and allows for unitary exchange of colors; 
it is described by the same unitary hopping matrix U x for both \gi) 
and |e;) states. Hopping along the y-direction is also laser assisted 
and attains "spin dependent" phase factors. Bottom) Trapping po- 
tential in y-direction. Adjacent sites are set off by an energy A due 
to the lattice acceleration, or a static inhomogeneous electric field. 
The lasers flu are resonant for transitions \gu) <-> \e2i), while Q^i 
are resonant for transitions between \eu) <-> \g2i) due to the offset 
of the lattice sites. Because of the spatial dependence of Q,i : 2 (run- 
ning waves in ztx direction) the atoms hopping around the plaquette 
get the unitary transformation U — Uy(m)U x U y (m + l)Ul, where 
U y (m) = exp(27rim diag[«i, 02]), as indicated in upper figure. 
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elements of our scheme, as shown in Fig. [2] are experimen- 
tally accessible. Nevertheless, consistent gauge group real- 
izations demand tunneling matrix amplitudes to be controlled 
in such a way that they strictly belong to the corresponding 
unitary groups, which assures the vector potentials to be com- 
posed of the Hermitian generatorsl 18]. The scheme of Fig. [2] 
allows to generalize the Hamiltonian Q to the case of non- 
Abelian vector potentials. In fact, we replace the components 
of A by the corresponding matrices from the group algebra. 
In particular, the illustrated setup generates "artificial" gauge 
potentials of the form A — (A x , A y (m), 0), with 
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More precisely, the scheme associates a unitary tunneling op- 
erator with every link in analogy with standard lattice gauge 
theory prescriptions Il2ll : U(m — l,n — > m, n) = U x , 
U(m,n — > m — l,n) = U\, U(rn,n — » m,n + 1) = 
U y (m), U(m,n — > m,n — 1) = U^(m), where U x = 
exp (+ieaA x /ch) and U y — exp (+ieaA y (m)/ch). The 
only difference is that A acquires an overall factor of hc/ea; in 



effect, though it does not behave well in the continuous limit 
a — > 0, the "magnetic flux" per plaquette, a± t 2 remains finite. 
Thus, we are in the same limit of ultra-intense fields ~ 1/a 2 
as in the "classic" Hofstadter case of Abelian magnetic fields. 
The ansatz ip(ma, na) — e lun g(m), leads to a generalized 
Harper wave equation 



g(m+l] 
g(m) 



= B{m) 



g(m) 
g(m-l) 



(5) 
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where e m (a, v) = e — 2cos(27rma — v) is the Harper en- 
ergy term. Both, ^(m, n) and g(m) are now two-component 
objects. In the particular case of eq. @, when two successive 
transfer matrices B(m) are multiplied by each other, eq. Q 
decomposes into a pair of equivalent 2D equations. Never- 
theless, to obtain the spectrum, one has to rely on numerical 
methods. Given each on — Pi/qi rational, the problem is Q- 
periodic (where Q equals the smallest common multiple of q\ 
and q2). The allowed energies are those for which the product 
of Q successive matrices B{m) has an eigenvalue of modulus 
1 . The spectrum shows a band structure, is continuous, and is 
bounded by two hyperplanes (depicted in grey in Fig. [3}- It 
exhibits a very complex formation of holes of finite measure 
and various sizes, which we name the Hofstadter "moth". Al- 
though a rigorous proof can not be provided, the "moth" re- 
minds of a fractal structure. Obviously, this fractal structure 
will be very sensitive to any sort of perturbation (finite size 
of the system, external trapping potential etc.) on very small 
scales. But, since the holes are true 3D objects with finite 
volume, the spectrum will be more robust on a larger scale to 
perturbations than in case of the Hofstadter "butterfly". This 
may be revealed in comparison of our "moth" with two uncou- 
pled Abelian "butterflies". The latter would form a spectrum 
of allowed energies consisting of intersecting lamellas of zero 
width. Here, the perturbations are required to guarantee al- 
lowed regions of finite volume. Since ol\ and a 2 are easily 
controllable, the spectral structure is directly measurable; a 
possible scheme would be to load a dilute (weakly interact- 
ing, or practically non-interacting) Bose condensate into the 
lattice and measure its excitations, e.g., by looking at the time 
evolution of the particle density, as suggested in Ref. Al- 
ternatively, one could load an ultra-cold polarized Fermi gas 
and measure its Fermi energy as a function of the number of 
particles. 

It is interesting to consider yet another effect that becomes 
particularly spectacular in the limit of ultra-intense gauge 
fields, and that can be measured in the proposed system: a 
non-Abelian Aharonov-Bohm effect, to be considered as an 
example of non-Abelian interference. In order to realize it, 
one should prepare, for instance, a weakly interacting Bose 
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FIG. 3: The Hofstadter "moth" spectrum. Eigenenergies e are plotted 
versus aj — Pi/qi,£ [0, = 0.5] (i = 1,2), where q, < 41 and 




condensate in a definite internal state V'o) around a location 
Pi . Then, the BEC (or parts of it) should be split by Raman 
scattering, and piecewise dragged (using, e.g, laser tweezers) 
to a meeting point P 2 on two distinct paths. These corre- 
spond to the unitary transporters U\ and respectively. A 
measurement of the density of atoms at P 2 will reveal a non- 
Abelian interferometer signal, i.e., it will detect the interfer- 
ence term ni nt oc ("fol^^i \ipo)- Choosing, e.g., the rect- 
angular loop, consisting of L X1 steps in the — x-direction, L 
steps in the ^-direction, L Xl +L X2 steps in ^-direction, L steps 
in — y directions, and finally L X2 steps in — ^-direction, we ob- 
tain for the gauge potential of eq. (0}: jii nt oc 1 if L Xl + L X2 
is even, and n mt oc (* | exp [±2maL(L xl + L X2 )]\ip ), if 
L Xl +L X2 is odd, where a = diag(ai, a 2 ). The signal is thus 
extremely sensitive to L Xl and L X2 . If one attempted to mea- 
sure the phase shifts by introducing obstacles on the y arms 
of this interferometer, the result would strongly depend on the 
x-coordinate of the obstacles, and the location of Pi, P2; a 
non-Abelian manifestation of the external gauge potential. 

Obviously, the properties of the considered system in the 
limit of ultra-intense fields are quite complex. Though, to 
get a better intuition concerning the scope of this scheme, 
it is useful to consider also the "continuum" limit a — > 
with Vo — ► a 2 /m. Then, the Hamiltonian becomes H = 
(p — e/cA) 2 /2m. A natural question to ask is what kinds 
of gauge fields of "normal", i.e., a-independent strength, can 
be realized? In other words: what are the possible "artificial" 
gauge potentials that can be created? In general, phase factors 
resulting from running wave vectors can be introduced for all 
tunneling matrices: A(r) = ^(M x + [M 2 (f )+M 3 (f )], Ni + 
[N 2 (j) + iV 3 (f )], 0), where Afj, N t are arbitrary (in general 
non-commuting), dimensionless, and a-independent matrices 
from, e.g., u(2), and I is the characteristic length on which A 
varies. Note that in order to realize non-Abelian interferom- 
etry, it is sufficient to fulfill [Mi, iVi] ^ and set all other 
matrices to zero. This may be achieved using simpler exper- 
imental arrangements as those discussed above, and will be 



the subject of further investigation. Furthermore, local dis- 
order may be introduced in a controlled way that allows for 
small fluctuations of the matrices Mj. In particular, disorder 
can be made annealed, i.e. it changes on a time scale com- 
parable with the relevant time scales of the system, and, thus, 
mimics thermal fluctuations. It can be of significant ampli- 
tude, provided it does not drive the assisting lasers system out 
of resonance. For example, wave vectors of running waves 
can strongly fluctuate and force the magnetic fluxes to fluc- 
tuate correspondingly. Even more complicated spatial depen- 
dences, e.g., piecewise linearity, of A are feasible by using 
static electric fields, laser induced potentials, etc. Additional 
lasers may introduce local, and in general time dependent uni- 
taries. Such transformations would generate arbitrary local 
temporal components of the gauge potential, A (x, y). Al- 
though for Yang-Mills fields in 2+ ID, this component may 
be gauged out adapting the Weyl, or strict temporal gauge 
11911 . the corresponding gauge transformations may introduce 
more complex spatial and temporal forms of the remaining 
two components of A. 

Finally, let us discuss, in which sense the proposed 2D 
scheme might be useful to study lattice gauge theories (LGT) 
in 2+1D. There, one uses the framework of Euclidean field 
theory, and various methods of statistical physics, such as 
Monte Carlo methods, to sum over all configurations of gauge 
and matter fields at finite or infinite temperatures. Thus, gauge 
fields are dynamical variables in LGT, but are obviously not 
in the proposed scheme. Moreover, our scheme is realized 
in real rather than in imaginary time, thus, only a limited set 
of configurations of gauge fields can be generated. Neverthe- 
less, the big advantage of our proposal is that given a gauge 
field configuration, the dynamics of matter fields in real time 
are given for free. When experimentally realized, the system 
may be used as a quantum simulator of matter fields in a given 
external gauge field. By generating other configurations, we 
may try to "mimic" the Monte Carlo sampling of LGT. Aver- 
aging over both, annealed disorder and quantum fluctuations 
should approximate the statistical average in LGT. However, 
this inevitably requires that generated configurations represent 
the characteristic or statistically relevant ones of correspond- 
ing LGT phases. 

Apparently, with the limited gauge field configurations ac- 
cessible in the proposed scheme, not too much can be done. 
Anyhow, at least some configurations may be generated that 
share characteristics of LGT phases, e.g., an area law ful- 
filled by Wilson loops in the confinement sector. In fact, the 
SU(2) gauge potential of eq. @ with fluctuating, but anti- 
correlated fluxes a\ = —a 2 yields an area law for Wilson 
loops in the xy-plane, provided the probability distribution of 
the fluxes is Lorentzian. However, it would be desirable to 
create configurations that exhibit other characteristics of the 
confinement phase such as appropriate distributions of cen- 
ter vortices, Abelian magnetic monopoles, instantons, merons, 
calorons etc. 

We expect that this program will lead to many fascinating 
results. It should be stressed that Yang-Mills theories in 2+1D 
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are in the center of interest in high energy physics, as they de- 
scribe the high temperature behavior of 4D models 12 ill . Re- 
cently, there has been progress in understanding these theories 
in the continuum limit, and in the pure gauge sector: 
the gauge-invariant degrees of freedom have been identified 
and the Hamiltonian has been constructed. The ground state 
wave function is, to a very good approximation, a Gaussian 
function of the currents, and quantities as the string tension 
are known exactly. Our lattice models may shed new light on 
these recent discoveries. 
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